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Abstract 



We bring forward rather simple algorithm allowing us to calculate the effective impedance of 
inhomogeneous metals in the frequency region where the local Leontovich (the impedance) boundary 
conditions are justified. The inhomogeneity is due to the properties of the metal or/and the surface 
roughness. Our results are nonperturbative ones with respect to the inhomogeneity amplitude. They 
are based on the recently obtained exact result for the effective impedance of inhomogeneous metals 
with flat surfaces [|l| |^ . One-dimension surfaces inhomogeneities are examined. Particular attention 
is paid to the influence of generated evanescent waves on the reflection characteristics. We show that 
if the surface roughness is rather strong, the element of the effective impedance tensor relating to 
the p- polarization state is much greater than the input local impedance. As examples, we calculate: 
i) the effective impedance for a flat surface with strongly nonhomogeneous periodic strip-like local 
impedance; ii) the effective impedance associated with one-dimensional lamellar grating. For the 
problem (i) we also present equations for the forth lines of the Pointing vector in the vicinity of the 
surface. 

' 

^ Keywords: Inhomogeneous metals, surface impedance 
O PACS number: 78. 90. -ft 

^ 1 INTRODUCTION 

>'- . 

QQ It is well known that to solve an electrodynamic problem external with respect to a metal, one can use 
04 ihe impedance boundary conditions jsj 

Ei = C[n,HJ, (1) 

Et and Ht are the tangential components of electric and magnetic vectors at the metallic surface; n is 
' ^ the unit external normal vector to the surface of the metal; ^ is the two-dimensional surface impedance 
tensor. Its real and imaginary parts define the dissipated energy and the phase shift of the reflected 
"j^ electromagnetic wave respectively (see, for example, Q). 

I If for an inhomogeneous metal the penetration depth of an electromagnetic wave S is much less than 
I a characteristic size a of the surface inhomogeneity, omitting the terms of the order of \zeta\5/a, the local 
Values of the impedance can be used in the boundary conditions (1). In this case Eqs.(l) are usually 
^ called the local Leontovich (the impedance) boundary conditions Q. The surface inhomogeneity can be 
Q due to the inhomogeneity of the properties of the metal or/and the surface roughness. The elements of 
• • the local surface impedance tensor of an inhomogeneous metal are functions of position on the surface. 
. ^ If the metal surface is rough, the components of the unit vector n depend on position too. The more 
complex is the geometry and composition of the surface, the more complex is the solution of the problem. 
5-H When the surface inhomogeneity is strong, the solution requires cumbersome numerical calculations. 
. I The structure of electromagnetic fields in the immediate vicinity of the inhomogeneous metallic surface 
is very complex too. The fields above the surface are represented by incident, reflected, scattered and 
evanescent waves. However, sometimes it is sufficient to know the reflected electromagnetic wave only, 
not going into details of the field distribution near the surface. 
In what follows, when calculating we assume that 

(5 < a < A, (2) 

where A = 2ttc/uj is the vacuum wavelength. When inequalities (2) are fulfilled, only the fields averaged 
over the surface inhomogeneities " survive" far (at the distance much greater than the vacuum wavelength) 
from a stochastically or periodically inhomogeneous metallic surface. To calculate these fields, it is 
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convenient to introduce the effective surface impedance relating the tangential averaged fields at a plane 
close to the real metallic surface (or coinciding with this surface, if the surface of the metal is flat). 

We define the effective impedance tensor with the aid of equations similar to Eqs.(l) for a homogeneous 
metal with the flat surface, let us say, X3 ~ 0: 

< Et >=Ce/[e3,< Ht >]; (3) 

63 is the unit normal vector to this plane, <:> denote an average over the plane x-i = 0. In the case 
of a stochastically inhomogeneous surface the averaging is carried out over the ensemble of realizations 
of the inhomogeneities; if the inhomogeneity is a periodic one, the brackets denote the average over the 
period. If the tensor C,ef is known, when calculating the averaged electromagnetic fields we return to the 
reflection problem from a homogeneous flat surface. 

The general way providing the correct result for the effective impedance is to derive equations for 
the averaged electromagnetic fields both in the metal and in the medium over the metal (we assume it 
to be vacuum), as well as the boundary conditions for these fields at the metal- vacuum interface. By 
proceeding from these equations the tensor can be calculated. 

The similar method is used rather often when calculating effective characteristics of inhomogeneous 
media. It goes back to the pioneering works of I.M.Lifshitz and his co-authorsEl However, mostly this 
regular method is applicable when the inhomogeneity amplitude is small and the perturbation theory can 
be used (with regard to the calculation of the effective impedance see [|[ 0, ^). The exact solutions for 
effective characteristics of strongly inhomogeneous media can be found very rarely. One of examples is 
the calculation of the effective static conductivity of some two- dimensional inhomogeneous metals |^ . 

The other example is a quite new result for the effective surface impedance of an inhomogeneous 
metal with flat surface |^. It is valid in the frequency region of the local impedance boundary 
conditions applicability. Recently this solution was used when calculating the effective impedance of 
strongly anisotropic polycrystalline metals both under the conditions of normal |^ and anomalous skin 
effect [|0|. 

In |l|r|§l the expression for the effective impedance was obtained with the aid of the general scheme. 
Herein we present much more straightforward reasoning leading to the same result and extend it to the 
case of rough metallic surfaces. This result is a nonperturbative one with respect to the inhomogeneity 
amplitude. It is exact up to the limits of the local impedance boundary conditions (1) applicability. 

The organization of the paper is as follows. In Section 2 some general remarks useful when calculating 
the effective impedance are given. In Section 3 the simple way to derive the expression for Ce/ in the 
case of fiat inhomogeneous surfaces is presented. As an example, we analyze in details the case of one- 
dimensional strongly nonhomogeneous periodic strip-like local impedance. In Section 4 we present the 
algorithm allowing us to calculate the effective impedance in the case of one-dimensional rough surface. 
The obtained formulas are used to calculate the effective impedance associated with a lamellar grating. 
The s- and p-polarization states are examined separately. Concluding remarks are given in Section 5. 



2 PRELIMINARY REMARKS 

Suppose the local impedance tensor C is a known. A plane wav SW = Eoe-'('=^='+'^*) is incident upon the 
surface from the region 13 > 0. The metal occupies the volume beneath the plane 0:3 = 0. The surface of 
the metal is assumed to be the plane X3 = itself, or it is the plane X3 = perturbed by a roughness. For 
the sake of definiteness we suppose that the surfaces does not stand out over the plane 0:3 ~ 0. Suppose 
the surface inhomogeneity is a stochastic or a periodic one. For a rough surface its profile is defined. 

Let E(x) and H(x) be the total electromagnetic vectors in the region X3 > 0. We extract the averaged 
electric and magnetic vectors: 

E(x) E(x) > +e(x); H(x) =< H(x) > +h(x); < e(x) >=< h(x) >= 0. (4.1); 

If a <C A, the fields e(x) and h(x) represent mainly evanescent waves damping at the distances of the order 
of a above the surface. When X3 > d, where d is a distance complying with the inequality a ^ d ^ A, 
the leading terms in the expressions for the total fields are equal to their averaged values: 

E(x) E(a;3) Eoe-**^"=^'+ < E'' > e'^^^; H w < H(a;3) >= Hoe-''=^^'+ < H'' > e''^^^; (4.2) 

^The first one was |Q 

^To calculate the surface impedance of a metal it is sufficient to investigate the normal incidence of an electromagnetic 
wave onto the metal surface (see, for example, 111)) 
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< E'' > and < H'' > are the electric and magnetic vectors of the reflected wave respectively. 

Suppose we are interested in the averaged fields only. Since no inhomogeneous parameters enter 
the Maxwell equations in the vacuum, the averaging of these equations provides the standard Maxwell 
equations for the averaged fields. Our problem is to obtain the boundary conditions for the averaged fields 
with the aid of the exact Eqs.(f). Comparing the last with Eqs.(3), we define the effective impedance. 
In what follows we solve this problem, reducing it to calculation of the magnetic field in the vicinity of 
infinitely conducting surface of the same geometry. 

To start the calculation, we remind that the elements of the tensor are of the order of the ratio 
S/X <C 1. Consequently, in accordance with Eqs.(l), the tangential components of electric vector at the 
metallic surface are much less then the tangential components of magnetic vector. 

Next, at the surface of a metal the magnetic vector Ht is approximately equal to the magnetic vector 
Hj'^'' at the surface of a perfect conductor {\Cap\ = 0) of the same geometry. The same is true in the 
immediate vicinity of the metallic surface. Then at the plane X3 = 

<H,>«<Hr(x||,0)>, (5) 

X|| is two-dimensional position vector at the plane 2:3 = 0. With regard to Eqs.(l) the first nonvanishing 
term in the expression for the tangential electric vector at the metallic surface is 

E,«C(x)[n(x),Hr]. (6.1) 

If the surface is flat (it is the plane X3 = itself; n = 63), the vector < Ej > is obtained from Eq.(6.1) 
directly: 

<Et> = <C(x||)[e3,Hr]>; (6.2) 

In we showed that the corrections to the expressions (6.2) are of the order of S^/uXEq. These small 
terms cannot be taken into account in the framework of the local Leontovich boundary conditions and, 
consequently, the averaged electric field cannot be calculated more accurately. 

In the case of a flat inhomogeneous metallic surface Eqs.(3), (5) and (6.2) allow us to calculate the 
effective impedance. The obtained result cannot be improved. It is exact up to the limits of local 
impedance boundary conditions applicability. 

The equations (5) and (6.1) allow us to define the effective impedance for a rough metallic surface 
too. However, the method of calculation is not so straightforward (see Section 4). 



3 FLAT SURFACE OF INHOMOGENEOUS METAL 

When the aforesaid argumentation is taken into account, the result of ||l|, |^ for the effective impedance 
of an inhomogeneous metal with a flat surface, 

Ce/=<C(X||)>, (7) 

is almost obvious. Indeed, in the case of normal incidence of an electromagnetic wave upon the flat 
surface X3 = of a perfect conductor the magnetic vector at the surface is H^'^'^ = 2Ho, where Hq is 
the amplitude of the magnetic field of the incident wave. This vector does not depend on position at the 
surface. Thus, with the aid of Eqs.(5) and (6.2) we have: 

< Ht 2Ho; < Et 2 < C > [63, Hq]. (8) 

Comparing Eqs.(8) with Eqs.(3), we obtain Eq.(7). This result is correct for any metal with inhomoge- 
neous flat surface both under the conditions of normal and anomalous skin effect. 



3.1 ONE-DIMENSIONAL STRIP-LIKE LOCAL IMPEDANCE 

Suppose the local impedance tensor of an isotropic inhomogeneous metal is one-dimensional periodic 
strip-like function: Qk{x\\) = Cixi)Sik (hk — 1,2), where 

C{xi) = Ci; |a;i|<a, 

C2; a < \xi\ < 6, ^ ' 
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2b is the period. All the parameters, namely, a, b, b — a and the local impedances is Ci, C2 are compatible 
with inequalities (2). 

We use this example not only to illustrate how simple the effective impedance can be calculated with 
the aid of Eq.(7), but also to visualize the difference between the fields in the immediate vicinity of 
one-dimensional inhomogeneous metallic surface for the s- and p- polarization states. 

With regard to Eq.(7) for both polarizations the effective impedance is 

Ce/=Cl^+C2(l-^). (10) 

However, the difference between the polarizations exhibits itself in the values of the fields near the 
surface. To complete the picture, for both polarizations we calculate the fields above the metal and the 
time- averaged Pointing vector S = c[E, H]/47r defining the time-averaged energy flux. The continuity 
equation for this vector, divS = 0, represents the energy conservation law (see, for example, [^). 
P-POLARIZATION 

Let the magnetic vector of the incident wave be (0, Hq, 0) exp[— ifcxa] . The only nonzero component 
of the total magnetic vector above the metal can be written as 

H2{xi,X3) = Ho {e-'"^' + hoe'''''' + h{xi,X3)} , (11.1) 

where ho is the amplitude of the reflected wave and 



h{xi,X3)= Yl V^"«"^/'-"'"^); ag^^J{q7T/bf-k^^^. (11.2) 



-oo;g^0 



represents the evanescent waves. With regard to the Maxwell equations the components of the electric 
vector are 

Ei{xuX3) = ^0 < ^ £ /i,a,e(*"?"i/''-"'"3) + /iqC^^'^"^) - e^-''^"^) I ; (12.1) 

00 

E3{xi,X3)^-Ho^^ gV^^"'"^^'^"'"'^- (12-2) 

— oo;q7^0 

We see that to calculate the components of the electric vector E(a;i, X3) up to the terms of the order of 
^ allowed in the framework of the local impedance boundary conditions applicability, we need to know 
the amplitudes of the evanescent waves hq [q ^ 0) up to the terms of the order {kb)C, ^ C- 

According to Eqs.(l) in the p-polarization state i?i(a;i,0) — ~C{xi)H2{xi,0). We Fourier analyze 
this equation making use of Eqs.(ll) and (12.1). Then, omitting the terms of the order of C^, we obtain 
the coefficients hq-. 

/io-l-2Ce/; hq\q^o = ^ika)iCi-C2)Fi7:qa/b), (13) 
where F(x) = sinx/x. As a result up to the terms of the order of ^ 



■.{x,y) = 2Ho {cosikby/TT) - Qje^^'^'y/^^] ; (14.1) 



Ei{x,y) = 2Ho |isin(fc5y/7r) - (efe'^'^'''^^^ - ^(Ci - C2)C^i(^, ; (14.2) 

Es{x,y) - ^(Ci - C2)5i(x,y). (14.3) 

TT 

We introduce x = irxi/b and y = nx^/h^ The functions Cm{x,y) and Sm{x,y) are respectively the real 
and the imaginary parts of the functional Rm of the complex argument z = x -\- iy: 

i?™W=f:^^^e-. (15) 
9=1 ^ 



■^It is easy to calculate the explicit form of the functions Rm.(z) for an arbitrary m, but trying to be short we do not 
present here the relevant expressions. 
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Note, all the corrections taking into account the finite conductivity of the metal are of the order of the 
input impedances Ci and C2- Next, the expression for the magnetic field, Eq.(14.1), is just the same as 
if calculated for the p-polarized wave incident upon a homogeneous metallic surface with the impedance 
Cef- However, the evanescent waves contribute to the electric field. In equations (14.2) and (14.3) they 
are represented by the terms including the functions Ci{x, y) and Si{x, y). These terms are exponentially 
small when y 1 (or 3:3 » 6). 

From Eqs.(14) it follows that for the p-polarization state up to the terms of the order of C, the Pointing 
vector S(J')(the superscript [p) indicates the polarization) can be written as 

S(p){x,y) = -^^^yFj,{x,y), (16) 

TT 

where 

FA^.y) = |ReCe/ - ~ cos{kby/7T)C2{x,y). (17) 

Far from the surface (y ^ 1) the component s[^^ of the Pointing vector is exponentially small, and the 
component sjf'' = — (c/27r)|i7o|^Rt;^e/- At the surface s[^^ = and the component S^^^ is defined by the 
local impedance: = -(c/27r)|i7o|^ReCi, if |xi| < a and S^^^^ = -(c/27r)|i?o|^R.eC2, if a < |xi| < b. 
It is easy to see that along with Eq.(16) the components of the vector S can be written as 

S<"^ = ^rotA,(a=,t/), (18.1) 

TT 

where the vector Ap = —(0, Ap, 0) and 

Ap{x, y) = ^ReCe/ + ^""^^^ ~ cosikbyMS2{x, y). (18.2) 

From Eq.(16) it follows that in the p-polarization state in addition to equation divS'J') = we also 
have rotS(p) = 0. Thus, in the p-polarization state the vector S(p) satisfies the same set of equations as 
the electric vector in two-dimensional electrostatic problems. Making use of this analogy we introduce a 
complex potential w{z) = Fp{x,y) — iAp{x,y) for the Pointing vector S^p): 

^(^) = I^ReCe/ + Re(Ci - C2)i?2(^)| . (19) 

At the plane {x,y) equations Rew(z) = const define the equipotcntial lines for the vector S(p), and 
equations Iinw{z) = const define the force lines of this vector. Evidently, far from the surface (y ^ 1) 
the force lines are parallel to the y axis. Near the surface {y < 1) they are distorted due to the influence 

of evanescent waves. 

As an illustration, let us examine the limiting case of very strong inhomogeneity: a/b <g; 1, IC1I/IC2I ^ 
1 and (a/6)|Ci| '> |C2|- With regard to Eq.(lO) the relevant effective impedance is Ce/ = {0'/b)(i, and 
according to Eq.(19) the complex potential is 

^(^) = ^^ReCef |-^^ + ln(l - e'^)| . (20.1) 

The force lines corresponding to this potential are given by equations 

tan(a;/2) = tan(xoo/2) " [ , (20.2) 
1 + e V 

where defines the position of the given lino far from the surface. These force lines are shown in Fig.l 
(full lines) . All the lines being distributed uniformly far from the surface, meet at the point a; = on the 
surface. Immediately near the surface the force lines are the straight lines y = ax with a — cot(a;oo/2). 

It can be easily verified that the complex potential defined by Eq.(20.1) is the same as the electrostatic 
complex potential of the system of charged straight lines that are perpendicular to the plane {xi,x^) and 
intersecting this plane at the points x\ = 2bn; n = 0, ±1, ±2 : .... 
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S-POLARIZATION 

Let the electric vector of the incident wave be = (0, -Eq, 0) exp[— ifcxa]. The total electric vector 
above the metal is 

oo 

E2{xi,xs)=Eo{e-'''^' +eoeikx3 + e{xi,X3)}, e{xuX3) = ^ ^^^i^^g^i/b-a,^,) ^ ^21) 

— oo;g7^0 

Uq is defined in Eq.(11.2); eo corresponds to the reflected wave and e{xi,X2) represents the evanescent 
waves. 

To define the Fourier coefficients eo and Cg we use the boundary conditions (1). Omitting the terms 
of the order of C^/kb, we obtain 

2a 

eo = -(1 - 2Ce/); e,|q^o = y(Ci-C2m7rga/^). (22) 

(compare with Eq.(13)). Then in terms of the variables x and y our result for the fields above the surface 
is 

E2{x,y) = 2Eo |-isin(A;6y/7r) + Q/e^'^y/^ + ^(Ci - C2)Ci(a;,y)| , (23.1) 
H,{x,y) = 2Eo |cos(fc6y/7r) - Qfe'^'y/^ - ^^%^<^o(x, y)| , (23.2) 

H3{x,y) = 4tEo^-^^^^So{x,y); (23.3) 

In contrast to the p-polarization state the evanescent waves contribute both to the electric and mag- 
netic fields near the metallic surface. Next, the contributions of evanescent waves to the components of 
the magnetic vector arc of the order of C,/kb ( (compare with Eqs.(14)). 

Of course, in accordance with the energy conservation law in the s- polarization state the components 
of the time-averaged Pointing vector satisfy the equation divS(*) = 0. However, the curl of this vector is 
not equal to zero and, consequently, in the s-polarization state no complex potential can be introduced. 

To define the equation for the force lines of the vector S^") we write it as 

I IT* 1^ 

SW = ^^rotA,(a;,y), A, = -(0,^,0); (24) 

TT 



A«(.,.) = fReCe/ + 5^%^ 



kb 

sm{kby/TT)Si{x,y) H cos{kby/TT)S2{x,y) 

TT 



(25) 



(compare with Eqs.(18)). Then the force lines of the vector S^*) are given by equations As = const. In 
the case of very strong inhomogeneity {a/b <^ 1, |Ci| ^ IC2I and (a/6)|Ci| '> \C2\), they are described by 
the equation 

arctanu — 2- — ^ ^ - — — — —— = — — u(x,y) = - — - — cot(a;/2) (26.1) 
(1 -h u2) (1 - e-2a) 2 ' ^ "^^ 1 + e-y ^ ' ' ^ ' 

(compare with Eq.(20.2)). Again Xrx. defines the position of the line far from the surface. Immediately 
near the surface the force lines are the straight lines y = uqx. For the given value x^o the slope uq is 
defined by the equation 



arctanuo — ^ = — . (26.2) 

1 + ui ^ 



Uq TT — X, 

These force lines are shown in Fig.l (dashed lines). 

Summarizing the results of this section we would like to emphasize that the main difference between 
the two polarizations is the presence of the terms of the order of C/kb ~ 6/b in the expressions for the 
components of the magnetic vector in the s-polarization state. 
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4 1-D ROUGH SURFACE 

In this Section we show that if the mctalhc surface is rough, to calculate the effective impedance it is 
sufficient to know the magnetic vector ilP<^^ near the surface of a perfect conductor of the same geometry. 
The main problem is to calculate the components of this vector. In the general case they depend on the 
surface geometry and, as a rule, on position at the surface. 

We restrict ourselves with one-dimensional rough homogeneous metallic surfaces only. To be definite, 
let us suppose that our one-dimensional surface is a periodic one, 26 is the period. A part of the surface 
is the plane ^3 = itself, and each period has a deepening with a throat of the length 2o (Fig. 2a). 

In the case of the p-polarization state we calculate the relevant element of the effective impedance 
tensor for an arbitrary surface profile. For the s-polarization state we perform the calculation for one- 
dimensional lamellar gratings. We use this example to show the difference between the polarizations. 

To show the method of calculation, let us suppose that we know the vector H^"^'' for a given surface 
profile. With regard to the definition (3) of Ce/ and Eq.(4.2) we have that far from the surface {x^ » A) 
the fields are the same as the ones above a flat metallic surface with the impedance equal to Ce/. This 
means that far from the surface the only nonzero component of the time- averaged Pointing vector is 

^3 = -^^^ReCe/ (X3»A), (27.1) 

where Hq is the amplitude of the magnetic field in the incident wave. On the other hand, since near a 
metallic surface the magnetic vector H is nearly the same as H^"^'', at the rough metallic surface itself 
the time-averaged Pointing vector is 



Stt 



Hr(x)| n(x), (27.2) 



n is the unit external normal vector to the surface (see, for example, js)). The energy conservation law 
requires divS = 0, and, consequently, the energy fluxes through the infinitely distant surface and the 
surface of the metal are equal. Then 

ReCe/ = -^ / |Hr(x)rd;, (27.3) 

the integration is carried out along one period of the line bounding the surface at the (xi, X3) plane. 

The last equation enables us to calculate the real part of the effective impedance. To calculate Im^e/ 
we take into account that though no absorption happens when an electromagnetic wave reflects from a 
rough surface of the perfect conductor, a phase shift takes place. In other words, some pure imaginary 
surface impedance can be associated with the rough surface of the perfect conductor. It is easy to 
understand that its leading term is of the order of kb (we remind that kb 1). Since \(^\ kS kb, in 
the case of the finite conductivity the same pure imaginary surface impedance defines the leading term 
in the expression for IvnC^ef- 

Now, supposing that C = 0, with the aid of Eq.(3) we can define ImCe/ as the ratio of the tangential 
electric and magnetic fields averaged across the plane 2:3 — 0. The result for the averaged tangential 
magnetic vector is rather obvious. Indeed, Eq.(3) have the form of the boundary conditions for a flat 
metallic surface. Consequently, the leading term in the expression for the averaged tangential magnetic 
field is < >= 2Ho. This statement is true both for a metallic surface and the surface of the perfect 
conductor. Next, when C = 0, only the throat of the deepening contributes to the averaged electric field: 

a 

< Ef >=YbJ Ef (a:i, X3 = 0). (28) 

— a 

If the magnetic vector HP'^'^ is known, the last integral can be calculated with the aid of the Stokes 
theorem or one of its modifications (see below) . 

Concluding this subsection we would like to note, that the components of the vector H^'^'' depend 
on the parameter fc6 -C 1 (see Eq.(2)). We show that when calculating the leading terms of ReCe/ and 
ImCe/7 we need to know the field Hp*^^ up to the terms independent of this small parameter only. 

In what follows p- and s- polarization states are examined separately. 
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4.1 CALCULATION OF (u^ (P-POLARIZATION STATE) 

Let us start with the p-polarization state (the magnetic vector H = (0,7?2,0) is parallel to the rulings). 
With regard to our notations the relevant element of the effective impedance tensor is Cn • 

It is easy to see that when a p-polarized wave is incident upon one- dimensional rough surface of 
a perfect conductor, and the characteristic sizes of the surface profile are compatible with inequalities 
(2), in the vicinity of the surface the expression for the magnetic field i?!""^ can be written as i?!"^'^ = 
2Ho + {kb)h{xi, X3; kb); kb <^ 1 . (In APPENDIX 1 we calculate H^'^'' for an infinitely conducting lamellar 
grating.) As far as the terms of the order of {kb)( falls outside the framework of the; local impedance 
boundary conditions applicability, the corrections {kb)h{xi,X3;kb) has to be omitted when calculating 
ReCe/ with the aid of Eq.(27.3). 

Next, to calculate the averaged tangential electric field < E^^'' >= (< e[^^^^ >, 0, 0) entering Eq.(28), 
we write the Stokes theorem: 

= (j) drW''. (29) 



y"dS[V,Ef 



The integral in the left-hand side of Eq.(29) is taken over the cross- section of the deepening (dS ~ 
(0,dS', 0)); the integral in the right-hand side is carried out in the clockwise direction over the contour 
and the throat of the deepening. Now, taking into account that on the surface of the perfect conductor 
the tangential electric field is equal to zero, with the aid of the Maxwell equations we obtain 



(xuxs). (30) 



The integral in the right-hand side of Eq.(30) is of the order of kSH2/b, S is the area of the deepening. 
Since according to Eqs.(2) the value of kS/b <C 1, when calculating the leading term in the expression 
for < (E^'^^^ > we again use the approximation H2'^'^{xi,X3) = 2Ho. 

Let I/O be the length of the surface contour with respect to one period. Now, for an arbitrary shape 
of the deepening we can write the expression for the effective impedance leaving only the leading terms 
of ReCe/ and Im^e/ : 

CS" = ,31) 

We see that independently of the shape and the size of the deepening Re^^/ > Re^ and lni<^i( , being 
of the order of k multiplied by a characteristic size of the surface profile, is much greater than Im^ ~ k6. 
From Eq.(31) it follows that in the p- polarization state an intensification of absorption due to the surface 
roughness, is merely geometrical effect relating to an increase of the area of the absorbing surface. 

Let us estimate the maximum value (max of ReCif • It is evident that (max corresponds to the 
maximum permissible value of the ratio Lc/2b. Taking account of the inequalities (2), suppose the length 
Lc is of the order of the vacuum wavelength A and the period 2b is of the order of the penetration depth 
S. Then, since |^| ^ S/X, we have Cmax ~ 1- (We would like to recall that for good metals |C| "C 1.) Of 
course, this is only the upper bound of the permissible value of Re^if • However, the estimate shows that 
when the length of the contour is much greater then the period (Lc ^ 26), the value of Re^ii exceeds 
ReC significantly. Such increase of the effective impedance has to manifest itself as giant absorption of 
the incident p-polarized wave. 

The surface with Lc » 2b can be realized, for example, in such a way. Suppose a planar metallic 
surface has periodic grooves. Suppose the boundary of the groove is a branching line (Fig.2.b). Let us 
say, this line has a fractal structure. Then the ratio Lc/2b can be very large and, consequently, the value 
of the ReC^j'''-' close to Cmax can be achieved. When the throats of the grooves are suflaciently narrow, 
almost all absorption takes place inside "the pockets" beneath the plane X3 = 0. 

4.2 CALCULATION OF cS^^ (S-POLARIZED WAVES) 

It is much more difficult to calculate the effective impedance relating to the s-polarization state (the 
electric vector E = (0, £"2,0) is parallel to the rulings). In our notations this is the element C22 of ^^le 
effective impedance tensor. The point is that for this polarization the magnetic vector H*"^'' varies signif- 
icantly in the vicinity of the surface. Inside the grooves the strength of the magnetic field exponentially 
decays when the distance from the plane xs = increases. Next, in contrast to the p- polarization 
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state, evanescent waves generated in the region xs > 0, contribute to the leading term of the tangential 
magnetic field Hf'^^. As a result even for rather simple surfaces, such as lamellar grating, the value of 

(e f) 

Q2 can be calculated only numerically. 

Let us write down the general expression for C22^'' ■ Suppose the surface of the metal is of the type 
shown in Fig. 2a, and we know the vector H^*^^ in the vicinity of the surface. Then we can calculate ReCe/ 
with the aid of Eq.(27.3). 

To calculate lm(ef we need to know the averaged tangential magnetic < Hf'^^ {xi,0) > and electric 
< £2"^^ {xi,0) > fields. Although the components of the magnetic vector depend on coordinates, with 
regard to the aforementioned arguments we can use the approximation < H^*^^ (xi, 0) >= 2iJo (Hq is the 
amplitude of the magnetic vector of the incident wave). To calculate < £2'^^ {xi,0) > we use the vector 
modification of the Stokes theorem for the vector E'"^'': 

y[[dS,V],E''^''] =_^[di-,EP^1- (32.1) 

The domains of integration are the same as when calculating the integral (29). Multiplying the vector 
equation (32.1) by the unit vector 83 with regard to the Maxwell equations and the boundary conditions 
for the vector EP^^ we obtain 

< En^uO) >= I ASHrixuXs), (32.2) 
Now it is clear that in the s-polarization state the the effective impedance is 

Cif=ReC^Z-'-^W, (33.1) 
where the leading terms of Z and W are 



1 



^ = TTTTT^ / |Hr(x)|'dl, W = — - / dSHAxuXs). (33.2) 
4:Lc\Hq\ J ZbUo J 

(compare with Eq.(31)). Again is the length of the contour bounding our surface with respect to 
one period, and S is the area of the deepening. The factors Z and W show the difference between the 
polarizations. 

4.3 CS^^ ASSOSIATED WITH LAMELLAR GRATING 

There are some regular methods allowing us to calculate the magnetic vector at the surface of a perfect 
conductor in the small roughness limit. However, if the surface roughness is not small, analytical solutions 
can be obtained for some specific surfaces only. Therefore, as an example, we examine an infinitely 
conducting lamellar grating (Fig. 2c). We assume that all the sizes of the contour, namely the period 26, 
the width of the rectangular grooves 2a, the depth of the grooves h as well as 6 — a, are in agreement 
with inequalities (2). 

The electromagnetic fields in the vicinity of infinitely conducting lamellar gratings have been examined 
by a lot of authors (see, for example, [^). Therefore we do not go into details of calculation, but 
present only the basic formulae used to obtain the result and some remarks relating to the computational 
procedure. 

To calculate the vector H^'^^, following |l^, we present the electric field in the central groove E2^^ = 
£^2 •'(xi, 2:3) (|a;i| < a; —h < 2:3 < 0) as series of modal functions (j)n\xi,X3,) that are the solutions of 
the Maxwell equations: 

00 

E^2'\xi,x:i) = ^B,^")0l")(xi,X3), (t>\t\xi,xs) = sin ^(xi - a) sin/3„(/i + 2:3), (34.1) 

n=0 

= fc^ — (7rn/2a)^ and Re/3„, Im/3„ > 0. The representation incorporates the boundary conditions for 
perfect conductors on the vertical and bottom horizontal facets of the groove. 



9 



In the region ^3 > we seek the field E^*^^ = E2^\xi, X3) as 

L q=~ca ) 

ao — k and aq — i\J (nq/b)'^ — fc^, if q ^ 0. The amplitude of the incident wave is equal to one, cq is the 
amplitude of the reflected wave and {q ^ 0) are assigned to evanescent waves. 

Taking account of the fields periodicity, matching the electric vector at the plane 0:3 = and the 
tangential magnetic vector at the throat of the central groove, after eliminating the amplitudes eg and e^, 
we obtain a matrix equation for the coefficients i32n-i- been repeatedly shown (see, for example, 

that in the case of the s-polarization state all the even coefficients i?2n are equal to zero. 

Let us introduce a dimensionless parameter 7 = fc5 <^ 1. From Eqs.(33) it follows that when calcu- 
lating Q2 we need to know the components of the vector 'HJ"^'^ up to the terms independent of 7 only. 
We seek the coefficients i32n-i ~ 1; 2...) as series expansions in powers of 7. The first non- vanishing 
terms of these series are proportional to 7. To calculate the magnetic vector up to the terms independent 

(s) 

of 7, we cut off the series at the first terms. In this approximation in place of i?2n-i introduce the 
coefficients F„: 

- (2n-l)sinhx„' ^" " 2a ' 
The coefficients Yn are the solution of the infinite set of equations 



E°° . . N coth Y„ 2 , , 



00 



A / N {qn)[l+ COS qnfi] 

A»™(m) =2^— 2... .2 ,^ -2T; ^^ = 2a/b. (36.2) 

q=i[m) - (2n - I) ][{qn) -(2m-l)] 

Next, with h['^\xi, X3) to denote H^'^^{xi, x^) in the region above the plane X3 = 0, and H*^~^(a:i, 0:3) 
to denote 'HP'^^{xi, x^) in the central groove, in terms of the coefficients Yn we have in the limit 7^0: 

ilj+'(a;i,0) = 2[l + (5iJ(a;i)], (37.1) 



-3=1 



Next, 



-.5-<--' -Tat^^^- (^(- - «)) (1 . f )) . (37.) 

n=l />- \ / 



Note, from Eqs.(37.1) and (37.2) it immediately follows that < SH{xi) >== and < H^{xi,0) >= 2. 

Thus, at first we need to calculate the coefficient Yn solving numerically the infinite set of equations 
(36). Then we use Eqs.(37) to calculate the components of the magnetic vector and define the real and 
imaginary parts of the effective impedance with the aid of Eqs.(33). We describe the computational 
algorithm in APPENDIX 2. 

Let us write the coefficient Z entering Eq.(33.1) as 

Z = Zpi + Zs, (38.1) 

where Zpi and Zg are related respectively to the parts of the plane X3 = incorporated in the contour 
and the inner facets of the grooves. With regard to Eqs.(37) we have 



^pi = JhTh)l \^ + SH{x^)\'dx^; 



(38.2) 



10 



^'^A{bThj\l '^^ ^(«'^3)fda:3 + i / \h[ \xi,-h)\\lx,y (38.3) 



Since the explicit forms of Zpi and Zs in terms of the coefficients Yn are very lengthy, they are not 
presented here. The coefficient W defines the imaginary part of the effective impedance. For our lamellar 
grating 

W=-^y ^" . (38.4) 

In Figs. 3 and Fig. 4 we display our numerical results for two different values of the grooves depth. We 
choose h/2a = 0.1 and h/2a = 5 to represent shallow and deep grooves respectively. The dots correspond 
to h/2a — 0.1 and the stars to h/2a = 5. 

In Fig. 3a we plot the ratio Re(2'2^^ /^^C versus a/b showing that the presence of the grooves leads to 
an increase of absorption. Our calculations show that if h/2a < 1, the value of M{a/b) increases with an 
increase of the ratio h/a. However, the results for M = M{a/h) are practically the same when h/2a > 1. 
This means that an incident wave "understands" the grooves, whose depth is more than 2a, as infinitely 
deep. In other words, for our lamellar grating the stars in Fig. 3a define the maximal values of the real 
part of the effective impedance in the s-polarization state. 

In Fig. 3b for the same values of h/2a we present the values of Z = 'Re(^22^'' /Ke(^[l'^'' for different ratios 
a/b. The coefficient Z defines the distinction between the values of the effective impedance for the s- and 
p-polarization states. In Fig. 3c we show the ratios Zpi/Zs versus a/b. This ratio defines the distribution 
of the absorbed energy between the horizontals sections 2:3 = and the inner surfaces of the grooves. In 
Fig. 4 the coefficients W — Im^^t'^'' /lni(^['j^'' versus a/b are plotted. 

Our results show that the same as for p-polarized waves, the presence of the grating leads to an 
increase of absorption of incident s-polarized waves. However, since Z < I, always the absorption of the 
p-polarized waves is more intensive. The more deep are the grooves, the more the difference between the 
polarizations. We also see (Fig. 4) that the same as the factor Z, the factor W entering the expression 
for the imaginary part of the effective impedance (see Eq.(33.1)) is less than one. However, in contrast 
to Z, it depends on the ratio a/b only slightly. 

We must note that a similar problem was solved in |0 by L.A.Vainshtein, S.M.Zhurav, A.I.Sukov. 
The authors of JT^ examined a one-dimensional grating composed of semi-infinite plates (infinitely deep 
grooves) exposed to s-polarized normally incident electromagnetic wave. They also made use of the 
impedance boundary conditions (1) and with the aid of some other method based on the solution of the 
Wiener- Hopf equation, calculated the reflection coefficient R. Then supposing that \C\/kb <C 1, in the 
framework of the perturbation theory the value of R was defined up to the terms of the order of C. The 
real part of the effective impedance was determined when calculating the difference 1 — \R\^ ■ 

Of course, two different methods of calculation must lead to the same results. However, the values 
of ReC22'^VR'6C for deep [h/2a > 1) grooves obtained in the present work (the stars in Fig. 3a) are 
approximately twenty percents less than the results of . Trying to find out the origin of this discrepancy 
and to verify our approach, in we examined the system investigated in using two different 
approaches. We reproduced the perturbation theory calculations of jl7| with the aid of the modal 
functions method used in the present work. The results were the same as in [ p^ . Then we repeated 
the calculation with the aid of Eq.(27.3). Comparing the results we obtained the same twenty percents 
difference. The analysis of showed that the only source of the discrepancy could be in application of 
the perturbation theory. 

Unfortunately, we were not successful in detecting accurately why the standard perturbation theory 
used in and |Q led to the results differing from the results obtained with the aid of our approach based 
on Eq.(27.3). However, we found some reasons allowing us to assume that in this problem the standard 
perturbation theory was not justified properly (see ||l^ ). We would like to note that the aforementioned 
difference is not of principle for qualitative description of the results obtained when analyzing the reffection 
of s-polarized waves from lamellar gratings. 

5 SUMMARY 

On the basis of the local impedance boundary conditions (1) we examined several metallic systems with 
strong surface inhomogeneities. For each system we calculated the effective surface impedance tensor. 
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Our results are exact within the accuracy of Eqs.(l). 

If the surface of an inhomogeneous conductor is flat, the effective impedance is equal to the values 
of the local impedances averaged over the surface (see Eq.(7)). As an example, we examined a flat 
metallic surface with one-dimensional periodic strip-like local surface impedance. In this case the effective 
impedance is an isotropic tensor: c\l^'' = CefSik («, fc = 1, 2); (ef is defined in Eq.(lO). 

We analyzed the role of evanescent waves paying special attention to the difference between the p- and 
s-polarization states. We showed that in the p-polarization state the evanescent waves did not contribute 
to the magnetic field in the immediate vicinity of the surface. On the contrary, in the s-polarization state 
just the evanescent waves provide the main contribution to the magnetic field taking account of the finite 
conductivity. 

For this problem we also calculated the components of the time- averaged Pointing vector S. For 
normally incident waves the force lines of the vector S being directed along the 0:3 axis far from the 
surface (x^ > A) become distorted near the surface. The stronger is the surface inhomogeneity, the 
more the distortion of the lines. We showed that in the p-polarization state the Pointing vector was 
defined by the complex potential w{z) (z ~ xi + 1x3) (see Eq.(19)). The aforementioned difference in the 
contribution of evanescent waves to the magnetic field is the reason why a complex potential cannot be 
introduced in the s-polarization state. 

When the inhomogeneity is very strong (narrow strips with rather large impedance separated by wide 
strips with very small impedance), for both polarizations all the force lines of the vector S meet at the 
surface in the narrow regions where the local impedance is big (see Fig.l). In the p-polarization state 
the equation for these lines is just the same as the equation for the force lines of the vector E in the 
electrostatic problem for the periodic system of charged straight lines. 

As an example of calculation of the effective impedance in the case of a conductor with rough surface 
we analyzed a periodic ID surface of homogeneous isotropic conductor depicted in Fig. 2a. In this case 
the effective impedance is a tensor whose elements differ significantly, when the grooves are rather deep. 

In the p-polarization state the electromagnetic field penetrates into the grooves on the surface. In our 
frequency region this means that the magnetic field £[2'^^ is nearly constant in the immediate vicinity of 
the surface. As a result, the surface as a whole takes part in absorption of p-polarized waves. This leads 
to a significant enlargement of absorption: the longer the length of the surface profile, the more the real 
part of the element of the effective impedance tensor relating to the p-polarization state (see Eq.(31)). 
In the framework of our approximation the upper bound for the real part of this element of the effective 
impedance tensor is of the order of one. The surfaces showing the giant absorption of p-polarized waves 
(Ce/ ^ 1) can be realized, for example, if the contour of the grooves is a strongly branching line. 

In the s-polarization state evanescent waves contribute to the magnetic field in the vicinity of the 
surface. As a result at the surface the components of the magnetic vector 'HJ"^'^ depend on coordinates. 
The most important that inside the grooves all the fields decay exponentially with the increasing distance 
from the plane x^ = 0. No simple expressions for this element of the effective impedance tensor can be 
obtained. For each surface the calculation has to be done separately beginning from the calculation of 
magnetic field near the infinitely conducting surface of the same geometry. 

As an example we examined the lamellar gratings shown in Fig. 2c. Recently many authors studied 
scattering from surfaces having rectangular grooves. These surfaces are of interest because they can be 
manufactured very easily providing a possibility to check theoretical results. Most of the investigations 
have been devoted to resonant enhancement processes (see, for example, |l3[ ^). However, the 
frequency region examined in these works does not comply with the inequalities (2) defining the framework 
of the effective impedance approach applicability. 

Our numerical results obtained when calculating the effective impedance associated with the lamellar 
grating for the s-polarization state are shown in Figs. 3 and Fig. 4. 

Although in the s-polarization state the presence of the grating leads to an increase of the real part 
of the effective impedance too, even for very deep grooves it cannot achieve such giant values that are 
possible in the p-polarization state. Taking into account that in the case of an arbitrary 1-D surface 
profile s-polarized waves do not penetrate into deep grooves, we can generalize our result: if an arbitrary 
1-D surface profile has rather deep grooves, the effective impedance associated with this surface is a 
strongly anisotropic tensor. Its element relating to the p-polarized state is much greater than the element 
relating to the s-polarized state. 
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6 APPENDIX 1 

Herein for the p-polarization state we calculate the leading term in the expression for the total magnetic 
field i?2*'" above the infinitely conducting lamellar grating depicted in Fig. 2c supposing that h,b <^ X. 

In the half-space above the grooves {xs > 0) we adopt the plane-wave representation of the magnetic 
field iff = ir2+(xi,X3): 

H+{xi,X3) ^ Hole-''''''' + f^+Qi{7.qx,/b+a,xs)\ ^^^^-^ 

[. q=-oo ) 

ttq is defined in Eq.(34.2); Hq is the amplitude of the incident wave. 

Inside the central groove (|a;i| < a; —h < x^ < 0) we seek i/f'^'^ = H2 {xi^x^) as a scries of modal 

functions (pn\xi,X3) that are the solutions of the Maxwell equations. With regard to the boundary 
conditions for perfect conductors we have 

00 

H^{xuX3) = HoY,B^f(t>lr\^uX3); 4^)(a;i,a;3) = cos(|^(a;i -a))cos/3„(/i + a;3), (A2) 

n=0 

Pn is defined in Eq.(34.1). 

Applying the boundary conditions on the plateaus X3 = and the continuity conditions across the 
central slit, we obtain the set of equations for the coefficients bI^^: 

00 

i sin(/3„/i)/(Pi - ^[S^) cos Pmh + (5„^oi?^^^ cos{kh)] = -20*^^, (A.3a) 

n=0 

00 ^ 



C7„, = 2qA{m; 2qa/b), A{m; z) = — ^ _ J2] ^' (^'^c) 

In terms of the coefEcients Bn^ the amphtudes h'^ are: 

00 

h+=Sgo + — J2 B^n^PnCnq sin /3„/l. (^.4) 

When the elements of the matrix /im are presented as the series expansions in powers of the small 
parameter 7 = kh, it is easy to see that the element Iq^ has the term a/267, and the expansions of all 
the other elements begin with the terms independent of 7. Then from Eqs.(A.3a) it follows that up to 
the terms independent of 7 

b'^p^ = 2; B^f) = 0, if n ^ 0. (^.5a) 
Within the same accuracy from Eq.(A.5) we have 

/ij = 1; h^+^ = 0, if g 7^ 0. (^.56) 

Thus, inside the grooves the modal function with n = gives the main contribution to the magnetic field. 

Next, the amplitudes of evanescent waves are negligibly small in comparison with the amplitude of the 
reflected wave. With regard to Eq.(A.l) and Eq.(A.2) this means that in the vicinity of the surface the 
magnetic field is nearly constant: 112^^ {x\,X3) w 2Hq. 
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7 APPENDIX 2 



In the s-polarization state when calculating the components of the magnetic vector H^'^'', special attention 
has to be drawn to the corner points of our lamellar grating. It is well known that if a surface has 
geometrical singularities, some components of electromagnetic vectors have singularities two. In the s- 
polarization state the tangential magnetic field in the immediate vicinity of a rectangular two-dimensional 
infinitely conducting wedge (the wedge is along the X2 direction) is proportional to p~^^^, where p — 
^Jx\+ x\ is the distance from the corner point (see, for example, ||l5|). Consequently, h'^\xi,Q) ~ 
Ixi ± ap"^^'^, when xi — > ±a, and H'^\±a,X3), when X3 — > 0. We need the solution of Eqs.(36) 
providing the aforementioned behavior of the magnetic field. 

In it was shown that the behavior of the fields near the edge is defined by the asymptotic behavior 
of the modal functions representing the fields. Let us show that in our case the tangential magnetic field 
near the edge is proportional to p^^^'^, if Yn ^ l/(2n — 1)^^"^, when n ^ 00. 

Indeed, let us examine, for example, \xitX^) when xi = a and X3 —0. It is the tangential 
magnetic field on the vertical facet of the groove near its throat. Taking into account that when X3 0, 
the behavior of h'^ \a,xz) is defined by the coefficients F„ with n ^ 00, with regard to Eq.(37.4) we 
obtain 

, 00 

ij(-)(a,X3)«^^r„e-"l-^l, 1^31-0. (A6a) 

n=l 

Suppose Yn ~ (2n — 1)'' when n — > cx). With regard to the asymptotic equality [p^ 

00 

^n^e""^ «r(l+p)z-(P+i), z ^ +Oand - Kp < 0, (^.66) 

n=l 

iT{x) is the Gamma function), comparing equations (A. 6a) and (A. 6b) we find that p = —2/3 has to be 
chosen. 

The next step is to search the solution of the set (36.1) truncating it to finite order N and increasing N 
until stability of successive solutions is obtained. Simultaneously, the sum defining the matrix A„m (see 
Eq.(36.2) has to be truncated to the finite order Q. The solution has to guarantee that Yn (2n — 1)"^^'^, 
when n ^ 1 . 

In [|l5j a similar set of equations was analyzed. This set was obtained when calculating the fields 
inside a rectangular one- dimensional branching wave-guide. It was shown that such set of equations 
allowed unrestricted number of solutions depending on the limit of the ratio Q/N for Q,N 00. The 
value of this limit ensuring the necessary behavior of the fields near the edge of the wedge determined the 
physically meaningful solution of the set in question. For a branching wave-guide this limit was calculated 
analytically as a function of the ratio of the sizes characterizing the wave-guide. 

Unfortunately, in the case of the lamellar grating we failed trying to define analytically the limit 
limAr,Q_^oo(Q/-^) providing Yn (2n — l)"^^'^ for n — > 00. Our results are based on numerical calcula- 
tions. 

In all the calculations we used N = 40. For this truncation number the stabilization was obtained 
in all examined examples. This value of N was sufficient to obtain a stable exponential solution Yn ~ 
((2n — 1) for 71 1 as well as the stable results for the real and imaginary parts of the effective 
impedance calculated with the aid of (33.1) and (38). 

To determine the ratio Q/N for given values of h/a and a/b, we solved the set of equations (36) using 
different truncation numbers Q and chose Q ensuring t = 2/3. We find that the ratio Q/N providing 
the given value of r, depends on the ratio h/a only slightly. When calculating Q/N as function of a/6 we 
saw that rather accurately the value r = 2/3 corresponded to the constant value of the ratio aQ/bN. For 
h/a — .2 and h/a = 2 the values of aQ/bN were 1.9 and 2.1 respectively. In Fig. 5 for a/b — l/sqrt{2) 
and h — 2a we show the exponent t as function of x: Q = [xbN/a] {[z] denotes the integral part of the 
number z). 

It worth to be mentioned that although for n ^ 1 the behavior of the coefficients Yn depends on the 
limAr,Q^oo(Q/-^)7 for small numbers n the coefficients Yn are practically the same when Qa/Nb > 1. 
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LIST OF FIGURES 

Fig.l. The force lines of the time-averaged Pointing vectors S(p) (fuU hnes) and S^") (dashed hnes) 
near the flat metaUic surface with strongly anisotropic strip-like local impedance. 

Fig. 2a. One-dimensional periodic surface that is examined in Section 4. 

Fig. 2b. A branching contour providing giant absorption of the incident p-polarized wave. 

Fig.2c. The lamellar grating configuration used to calculate the effective impedance in the s-polarization 
state. 

Fig. 3a. The ratio ReC22'^V-R'6C versus a/b for h = 10a (the stars) and h = 0.2a (the dots). 

Fig.3b. The values of Z versus a/h for h = 10a (the stars) and h = .2a (the dots). 

Fig.S.c. The ratio Z^ijZg versus ajh for h = 10a (the stars) and h = .2a (the dots). 

Fig. 4. The function W versus a/b for h = 10a (the stars) and h = .2a (the dots). 

Fig. 5 The exponent r as function of Q/N of x fov h = 2a and a/b = l/sqrt{2) {Q = [xbN/a]). 
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